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Abstract. A simplified model for the stationary, axisym- 
metric structure of magnetized winds with a polytropic 
equation of state is presented. The shape of the mag- 
netic surfaces is assumed to be known (conical in this pa- 
per) within the fast magnetosonic surface. The model is 
non-self-similar. Rather than solving the equilibrium per- 
pendicular to the flux surfaces everywhere, solutions are 
found at the Alfven surface where it takes the form of 
the Alfven regularity condition and at the base of the 
flow. This constraints the Transfield equilibrium in that 
the Alfven regularity condition is imposed and the regu- 
larity of the magnetic surfaces at the Alfven critical sur- 
face is ensured. The model imposes criticality conditions 
at the slow and fast magnetosonic critical points using 
the Bernoulli equation. These Alfven regularity and crit- 
icality conditions are used to evaluate three constants of 
motion, the total energy, angular momentum, and the ra- 
tio of mass to magnetic flux a, as well as the shape of 
the critical surfaces. The rotation rate f2 and the poly- 
tropic constant Q as a function of the magnetic surfaces, 
together with the mass-to-magnetic flux ratio on the axis 
ao entirely specify the model. Analytic results are given 
for limiting cases, and parameter studies are performed 
by numerical means. The model accepts any boundary 
conditions. Numerical calculations yield the value of the 
rotation parameter lo. Rotators can be deflned as slow, in- 
termediate or fast according to whether uj is much less or 
close to unity or near its maximum value for fast rotators, 
(|)2. Given the properties of astrophysical objects with 
outflows, the model allows their classification in terms of 
the rotation parameter. Critical surfaces are nearly spher- 
ical for slow rotators, but become strongly distorted for 
rapid rotators. The fast point remains at a finite distance 
for finite entropy fiows, in contrast to cold fiows.It is found 
that for a given mass loss rate, the rotation rate is limited. 
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1. Introduction 

Outflows are observed in a large class of astrophysical 
objects, from winds emanating from stars of all spectral 
types, to well collimated jets originating from young stars, 
compact objects and active galactic nuclei. Magnetically 
driven models are among those considered most promising 
for describing the generation of jets. In this case a large 
fraction of the stellar or accretion disc's rotational energy 
is converted into electromagnetic energy, and subsequently 
into kinetic energy. The influence of the magnetic field, 
and of the rotation on the flow, has to be assessed even 
for other types of objects in which the magnetic field plays 
a secondary role in the acceleration of the flow. The aim 
of this paper is to analyze the structure of outflows from 
rotating magnetized objects. 

Magnetized outflows start close to the central object. 
The mass loss rate is determined by the magnetic field 
configuration in this region and by the boundary condi- 
tions. After a major acceleration phase, the coUimation 
phase roughly starts at about the Alfven surface where 
the flow is deflected towards the axis by the so-called hoop 
stress. The problem of determining the stationary two- 
dimensional structure of the coUimation region of mag- 
netohydrodynamical outflows has not been solved. It re- 
quires the solution of the equilibrium of forces perpendicu- 
lar and parallel to the magnetic surfaces. One can describe 
the former by using the transfield or Grad-Shafranov equa- 
tion and the later using the Bernoulli equation for a poly- 
tropic equation of state. This situation is complicated by 
the existence of three critical surfaces which determine 
three constants of motion. Their position is only obtained 
as part of the global solution. 

Various approaches have been applied to address this 
problem. One approach assumes the shape of the field lines 
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to be given, thus ignoring the transverse equihbriuni (We- 
ber fc Davis ( |1967| ), Mestel ( |1968| ) Belc her & Me Gregor 
( |1976D Har tmann & Mc Greg or (|1982D , Pudritz fc No r- 
man ( [l98^ ), Mestel & Spruit ( |1987D , CassineUi ( [19891 )). 
The Bernoulli equation is then solved, together with two 
constants of motion, giving the mass to magnetic flux 
ratio and the total energy, which are obtained from the 
slow and fast magnetosonic critical surfaces. This allows 
the asymptotic speeds of the outflow to be determined, 
but it does not give any information on the coUimation. 
A variation of this approach are perturbative treatments 
of a spherically symmetric sol ution used for numerical 
simulations (Suess & Nerney ( 1973| )). For example, the 
treatment of an extended version of the Weber-Davies 
model in the space surrounding an axisymmetric system 
Sakurai (p^85| ),( |l987| ) and Uchida fc Shibata ( |l985| ) used 
a fully iterative numerical method to solve the Bernoulli 
equation together with the transversal force balance. 

Another approach uses the assumption of self- 
similarity where some specific dependence of the flow vari- 
ables on the independent variables is assumed (Chan fc 
Henrikse n (|1980D , Blandford fc Pay ne (|1982[ ), Pell etier fc 
Pudritz ( |1992| ), Tsi ngano s fc Sauty ( |1992D , S auty ( |1993D , 
Sauty & Tsinganos ( |1994D, G uyed fc Pudritz ( |1997| ), Gon- 
topoulos fc Lovelace ( |1994 ), Henriksen fc Valls-Gabaud 
(1994), Tsinganos et Trussoni ( |1991| ), Fiege & Henriksen 
(1996)). These models account for the force balance, but 
usually are not regular or valid in all space, and do not 
properly account for the fast magnetosonic critical sur- 
face. The axis of symmetry often appears as a singularity 



A more general approach is concerned with rigorous 
theorems on the asymptotic structure of magnetized winds 
(Heyvaerts fc Norman (1989)), where it has been shown 



for the electrical current (Blandford & Payne (1982), Pel- 
letier fc Pudritz ( |1992| )), though Pelletier fc Pudritz found 
a solution where the poloidal current did not diverge along 
the pole and at infinity. These solutions correspond to a 
current free plasma, where all the necessary poloidal cur- 
rent is concentrated along the polar axis. One should also 
mention the va riational approach presented by Rosso fc 
Pelletier (1994), and the slend er jet approximation devel- 
oped by Koupelis fc Van Horn ( |1989| ) and Koupelis ( |1990D . 



The disk-wind connection has also been central in 
many explanations of the origin of outflows. In some the 
engine responsible for the emission of the wind is a keple- 
rian disk threaded by a magnetic field that is either gener- 
ated in situ or advected-in from larger scales (Blandford fc 
Payne ( |1982| ), Kon igl (|1989| ), Pellet ier fc Pudritz ( |1992D , 
Wardl e fc K onigl ( |1993[), Li ( |1995D , Ferreira & Pelletier 
( |1993D ,( |1993| ), Ferreira ( [l99"7| )). In another type of model, 
a X-wind is postulated; the interaction of a protostar's 
magnetosphere with its surrounding disk opens some of 
the magnetospheric field lines to create the magnetized 
steU ar win d (Shu et al.( |198^ ), Shu et al. ( |1994| ), Najita fc 
Shu (1994)). Another method consists in studying numer- 



ically time dependent evolution o f the interaction between 
t he so urce and the flow (Guyed ( 1995 ), Guyed & Pudritz 



that outflows either become cylindrical at large distances 
from the source or parabolic, depending on whether they 
carry an electric current to infinity or not. 

This is intended to be the first of a series of papers 
on the structure of MHD outflows. We propose a simpli- 
fied model based on the assumption that the magnetic 
surfaces possess a shape in the collimation zone which is 
known a priori. But unlike the Weber-Davis type models, 
the balance of forces perpendicular to the magnetic sur- 
faces is taken into account on the Alfvcn surface, where it 
takes the form of the Alfven regularity condition, and at 
the base of the fiow. This constraints the transfield equi- 
librium in that the Alfven regularity condition is imposed 
and the regularity of the magnetic surfaces at the Alfven 
critical surface is ensured. Gnce given two constants of the 
motion describing the rotation and the thermodynamics 
close to the source as boundary conditions, the system of 
equations allows the determination of the three last con- 
stants of motion for otherwise general conditions. In the 
second paper, we examine coUimated outflows using these 
integrals of motion to determine the asymptotic structure 
of the flows. The third paper will address the question of 
jet stability with respect to magnetic instabilities which 
are formed in this model. 

This particular paper has a component devoted to 
our model and its astrophysical applications, and a sec- 
ond part which examines this model in context of other 
studies of outflow from magnetic rotators. The first part 
formulates the problem (§2), the assumptions and discuss 
the governing equations and relevant boundary conditions. 
We derive the new set of equations describing our model 
and obtain their solution analytically and numerically for 
slow (§3), fast (§4) and intermediate rotators (§5) and sug- 
gest astrophysical applications. The second part provides 
a comparison between studies and assesses the ability of 
this model to allow the classification of rotators (§6) in 
terms of suitable dimensionless parameters. It concludes 
with a summary which highlights the inner structure of 
magnetic rotator outflows (§7). 

2. Properties of MHD winds 

This paper deals with stationary and axisymmetric mag- 
netized rotating winds described in the framework of ideal 
MHD. 



Basic relations and notations Throughout the paper, 
{r,9,z) denotes cylindrical coordinates around the rotation 
axis while R stands for the spherical distance centered on 
the wind source. The vectors Br, eg, e^ are the orthogonal 
unit vectors associated with these coordinates. It is con- 
venient to split the magnetic field and the velocity into a 
poloidal part, which is in the meridional (r,z) plane, and 
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a toroidal part. The former is denoted by a subscript p 
while the latter is just the azimuthal component, so that 



B = B, 



Beeg 



(1) 



The poloidal part can be expressed in terms of a flux func- 
tion a{r, z) proportional to the magnetic flux through a 
circle centered on the axis passing at point r, z, as: 



1 9a 1 9a 

r oz r or 



Bb = T^Gr + --^€ 



(2) 



From this it follows that field lines of the poloidal field 
are lines of constant a, and so are the magnetic surfaces 
which are the surfaces generated by rotating them about 
the polar axis. 

The acceleration of gravity, g derives from a gravita- 
tional potential G{r, z) by g = — V G. The fluid velocity 
is denoted by v. We assume the density p to be related 
to the pressure P by a polytropic equation of state. This 
assumption replaces consideration of energy balance and 
is meant to represent simply some more complex heating 
and cooling processes. Then, we have: 



P = Qp-' 



(3) 



where Q is a factor related to the entropy of the flow and 
7 is the polytropic index. Q is constant following the fluid 
motion, but it can vary from one flow line to the next. 
From the stationary induction equation it results that flow 
lines follow magnetic surfaces. Therefore Q is constant on 
them and is therefore a function of a. 

A number of equations of stationary axisymmetric 
ideal MHD can be integrated to a set of equations ex- 
pressing the conservation of first integrals following the 
motion, namely the magnetic surface rotation rate O, the 
mass flux to magnetic flux ratio a on a magnetic surface, 
the specific energy E and the specific angular momen- 
tum of escaping matter, L. By the isorotation law and the 
specific angular momentum conservation law, the toroidal 
components of the velocity and of the magnetic field can 
be expressed in terms of the density and radius as 
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(5) 



These quantities would be singular when the denomina- 
tors vanish, when p = ^j-qO^, unless r^ becomes equal to 
^ when this happens. It can easily be checked that the 
poloidal velocity becomes in this case equal to the local 
Alfven velocity, calculated with the poloidal component 
of the magnetic field, vpA- For this reason /ioa^ is named 
the Alfven density 



Pa = Mo" 



(6) 



and L/il is the square of the so-called cylindrical Alfven 
radius at this Alfven point, 

r\ . ^. (7) 

More generally the subscript A will refer to values at the 
Alfven point and it can be shown that the alfvenic Mach 
number. Ma = vp/vpA, is given by 



^ P 



(8) 



The Alfvenic poloidal speed at the Alfven point can be 
expressed as 



VpA 



a|Va| 
pr 



(9) 



Basic equations The projection of the equation of motion 
on Bp yields by integration an equation, the Bernoulli 
equation, 

(10) 
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This equation can be also expressed, using equations (H) 
and (||) as 
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r'A PA- p : 

In the absence of MHD forces, this first integral would 
express the well known Bernoulli's theorem, i.e. the con- 
stancy of the sum of the kinetic, enthalpy and gravita- 
tional energy fluxes. The presence of the magnetic field in- 
troduces another energy flux, the Poynting flux, the fourth 
term in the equation. 

We introduce the following dimensionless quantities 
that depend on the particular ficldline defined by the mag- 
netic flux a: 



— The rotation parameter 
£trA 



ijj{a) = 



VPA 



(12) 



— The thermal parameter, closely related to the usual (3 
parameter at the Alfven point 



Pia) 



27 QpV 



(7 - 1) vpa'^ 
— The wind energy parameter 
2E 



6(a) ^ 



VpA-" 

— The gravity parameter 
2GM 



5(a) 



Ravpa 



(13) 



(14) 



(15) 
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These parameters weight the importance of associated ef- 
fects in terms of the magnetic energy at the Alfven point. 
/? is a measure of the importance of thermal effects in 
driving the outflow. 

The rotation parameter wiU play a crucial role in this 
problem and need to be clearly defined and explained. 
First one should take care about the difference between 
the definition given above and a different one used for ex- 
ample by Sakurai ( |1985| , |1987| ) or Spruit ( |1994| ) which is 



given by a; = 



^rA 



= . The limiting value 



3/2 



for 



^/GM/rA 

corresponds in both cases to the fast rotator limiting case. 
It is the minimum energy solution found by Michel (1969) 
for the classical approach of the Weber & Davics (1967) 
problem and gives a fast critical point rejected to infinity. 
With the latter definition omega will always be larger than 
this limiting value. The range of variations of co is rather 
different with our definition and need some explanations. 
For fast magnetic rotators, it can be easily shown in the 
frame of the Weber & Davies (1967) study that, neglecting 
the thermal effects, i.e. Q ~ 0, with respect to rotational 
and magnetic effects, two solutions exist only for uj larger 
or equal to the fast rotator limiting value (|) .If the 
thermal component of the flow is not neglected, the ro- 
tation parameter to still has the same limiting value but 
can become smaller than (|) and even reach zero in 
the case of vanishing rotation. This agrees with the re- 
sults given by Ferreira ( |1997D . Moreover in our model the 
rotation parameter lu increases for a given non-vanishing 
entropy as a function of the angular velocity and still has 

the same limiting value of (|) corresponding to the 
very fast rotator limit. There is actually no contradiction 
between these results and the cold wind theory since our 
limiting value is compatible with the inequality imposed 
by = 0. 

The projection of the stationary axisymmetric equa- 
tion of motion perpendicular to the magnetic surfaces ex- 
presses cross-field force balance in meridional planes. It is 
called the Grad-Schliiter-Shafranov equation or transfield 
equation (Okamoto ( |1978[ ), Heyvaerts & Norman ( |1989D ) 
and can be written as: 
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this point is the so-called the Alfven regularity condition, 
which is the particular form assumed by the transfield 
equation at the Alfven point. It can be written as (see 
Heyvaerts & Norman (1989) for details) 
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In this equation, p is the slope of the solution of the 
Bernoulli equation at the Alfven point, hereafter called 
the Alfven slope 



p=l 



dlogp 



^ d log r ^ 
that can be expressed as 

p=l 
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(18) 



(19) 



MHD flows have two other critical points which are 
brought about by the Bernoulli equation, the slow and 
fast magneto-sonic points which are located where the 
poloidal velocity equals one of the two magneto-sonic 
mode speeds. Any regular solution of the Bernoulli equa- 
tion must pass these points. Quantities referring to the 
slow or fast magneto-sonic critical point will be indicated 
by subscripts s or f respectively. Let the Bernoulli func- 
tion B{r,p) be the function defined by Eq. (fiVj), that has 
a constant value E{a) following the flow on this magnetic 
surface. The slow and fast critical points on a given mag- 
netic surface of flux parameter a are located where the 
differential of B{r, p) vanishes. The vanishing of the differ- 
ential form of the Bernoulli equation at constant a with 
respect to p and r stands as the criticality conditions 
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(21) 



rj , ^ So the set of equations is composed of the Bernoulli equa- 

tion and the Alfven regularity condition together with four 

or r c*^ (<]-iticality conditions, two deflned at the slow point and 
two at the fast magnetosonic point. 
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1. A model for the structure of rotating MHD winds 



It determines the shape of magnetic surfaces. Primes de- 
note derivatives with respect to a, i.e. E' = dE/da. This 
is a quasi-linear partial differential equation for the mag- 
netic flux function a{r,z). Eq . (p6 ) becomes critical at 
the Alfven surface (Sakurai (1985), Heyvaerts & Nor- 



man ( 1989 )) where it looses all its highest order deriva- 
tive terms. The condition for its solution to be regular at 



The main idea of our simplified model is to sacrifice exact- 
ness for simplicity, while retaining a high degree of gener- 
ality. Our model does not impose cross-field force balance 
everywhere but only at a few important places. First at the 
Alfven surface, where it takes the form of the Alfven reg- 
ularity condition. Second at the basis of the flow, which in 
the case of a low gas to magnetic pressure plasma amounts 
to assuming a uniform flux distribution on a small sphere 
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surrounding the wind source. Finally transfield force bal- 
ance can be imposed at infinity. This last point will be 
done in the second paper, allowing the problem to be more 
self-consistent. We solve the criticality conditions at the 
slow and fast mode critical points. Three first integrals of 
the motion are then determined from these conditions, and 
the Alfven regularity condition. The boundary conditions 
determine the other two. 

We moreover simplify by assuming the shape of 
poloidal field lines to be known up to the fast magne- 
tosonic point. We postulate, as a first approximation and 
for analytical convenience, that the deviation from con- 
ical shape is small enough in the region closer to the.JS 
wind source than the fast critical point to consider that 
the three critical points are aligned on conical magnetic 
surfaces. The shape of magnetic surfaces can indeed be 
regarded as conical at distances much smaller than the 
Alfven radius since the kinetic energy of the wind is in- 
sufficient to distort the magnetic field. In reality the mag- 
netic and kinetic energy start to become comparable at 
the Alfven point and the flow shapes the field at the fast 
point, especially if the fast point happens to be distant 
from the Alfven point. We do not assume, though, that the 
magnetic surfaces remain conical beyond the fast critical 
point. It is of course a first order approximation as can be 
seen in Sakurai ( 1987| ) where realistic magnetic field lines 
are obtained self-consistently by iterative method starting 
from a conical geometry. However this strong assumption 
will be relaxed in future works where more general forms 
of magnetic surfaces will be considered. One should con- 
sider this approximation as a first step in the study of the 
problem. 

Consequently a given magnetic surface is represented 
in this region by the equation z = r tan 0(a) and then 
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For a uniform distribution of the flux, the relation between 
the angle of a magnetic surface to the equator and the flux 
function a becomes 



cos^Q 
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(23) 



Then one can derive all the quantities referring to a mag- 
netic surface a defined in Eq. (p[),(p^), ( p^ and (p^), only 
in terms of the functions fi, Q, a, and of the distance to 
the origin and of the density at the Alfven point, Ra and 
PA- 
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7-1 



Fig. 1. A schematic representation of the structure of the 
model showing the conical shape of magnetic surfaces in 
the poloidal plane as light solid lines. The locus of each of 
the three critical points defines the corresponding critical 
surfaces. S, A and F refer to the slow magneto-sonic, the 
alfven and the fast magneto-sonic critical surfaces respec- 
tively. 
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Our set of equations consists of seven equations namely 
the Alfven regularity condition, the four criticality condi- 
tions defined at the slow and fast magneto-sonic points, 
and the Bernoulli equation written at the fast and slow 
critical points. The model is completely determined by 
the specification of the functions il(a), Q{a) and a on the 
axis (ao). The latter variable ap is directly related and al- 
most proportional to the total mass loss rate, and Q and 
Q, are given by the boundary conditions. 

The seven variables which appear in the system of 
equations are Rs^aj-, Ps,A,f and E. The positions and den- 
sities at the two critical magnetosonic points i?s, Rf, ps, 
and pf are given by the four criticality equations. The en- 
ergy E can be eliminated between the two Bernoulli equa- 
tions at the fast and slow point. The resulting equation 
and the Alfven regularity equation determine together pA 
and Ra- The energy E can be calculated in terms of quan- 
tities defined at the critical points thanks to the remaining 
Bernoulli equation which takes the following form 



A2 



(26) E 
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where the Re and pe are the positions and densities given 
at the fast or slow critical points. The positions and den- 
sities at the Alfven point are equivalent to L and a via 
Eq.(p|) and Eq.(ra) and can be associated to the last two 
equations of the system. Thus the three first integrals E, 
a and L are obtained. 



2.2. Analytical method 

2.2.1. General procedure 

We can obtain in important limit cases analytical expres- 
sions for the dimensionless specific energy, e and for the 
gravity parameter g as functions of the parameters il, Q 
and ao- We deduce the expressions for the positions and 
densities of critical points. As e and g have definite values 
on a given field line, the expression of these quantities as 
obtained from the slow point criticality relations and from 
the fast point criticality relations must be the same. These 
constants arc functions of i?s, Rf, ps, pf and, of course, 
of L and a on magnetic surface a. We use the following 
dimensionless variables 

Xg — ^ Xf — 

TA TA 



2 



Ps Pf 

Vs^ — , yf = — ■ 

PA PA 



(30) 



In the set of equations given in the previous part and defin- 
ing the model, one can calculate the Bernoulli equation at 
the critical points. First at the slow point. From Eq.([lC)|) 
this gives: 



e = Pvl-'+^H'i-ixl 
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Then at the fast point (in this case index s is replaced by 
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The positions and densities at the critical points are 
obtained from the vanishing of the differential of the 
Bernoulli equation (cf. Eq.^ and Eq.|l|). This gives the 
following four equations 
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(36) 



For the isothermal case, the equations (|33|) , (|34D for g re- 
mains the same and only the first term on the right hand 
side of the equations (|3l|),(p2|) for e are different. f3y]~^ 
and PyJ^ become respectively /31nj/s and f3hiyf. And 
the last two criticality equations are given by 
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(29) 2.2.2. A general constraint 



A condition that every solution has to fulfill can be derived 
from Eq.(pq), where all the terms must be positive. Since 
(1 — y/) is positive whereas (1 — j/s) is negative this implies 
that 



^(7-l)xlyri>l>^(7-l)4yJ+\ 
From this we derive the constraint 
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(For the isothermal case the constraint is given by r]:Pg > 
r'ip^r). This inequality can also be written in terms of phys- 
ically significant quantities as 



n-'p. 



PA 



(32) that is 



[S^H\ 



slow 



> 



[S^H] 



-Qp 



:7-l 



> 



n^Pfrj 



PA 



jQpr 



.(41) 



fast 



(42) 



where S is the Poynting specific energy and H the en- 
thalpy. Conditions arc thus found on the energies and on 
the positions and densities of the fast and slow magnetic 
point. The product of the square of the Poynting flux en- 
ergy by the enthalpy at the slow point has to be larger 
than the one calculated at the fast point. On the other 
hand, given the position and the density at the slow point, 
the possible values at the fast point are restricted by this 
inequality. 
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2.3. Numerical method 

Numerical investigations can study all ranges of parame- 
ters. The variables selected to describe the system remain 
r^, rf, TA, and Ps, Pf and pA with two free functions, Q(a), 
ft{a) and a free parameter, ao, the value taken by a{a) on 
the polar axis. The different forms of Q{a) and i^{a) that 
have been considered are given in Appendix B. It is con- 
venient to convert the algebraic magnetosonic criticality 
conditions into differential equations. The critical points 
can be followed from one magnetic surface to the next by 
differentiating these criticality equations with respect to 



d fdB 



d fdB 



da \ dr 



da\dr J ,■ 



da \ dp 



da\dp ) f 



(43) 



(44) 



Two more equations close the set, namely the Alfven regu- 
larity condition Eq.(n7h and an equation which expresses 
the fact that the specific energy has a given value on a 
magnetic surface, so that E{a) = B{rs,ps) — B{rf,pf) 
which can be expressed in differential form as 



-{B{r,,p,)-B{rj,pj)) 



0. 



(45) 



The system then consists of six differential equations for 
six variables. In practice, equations (Q) have been mul- 
tiplied by a factor (1 — — ) to prevent singularities near 
the axis where pf is almost equal to pA, since this fac- 
tor appears in the denominators. A similar problem arises 
close to the axis and on the axis in the calculation of the 
Alfvenic slope defined by Eq.dlS). We use the definitions 
given by Eq. (m) and Eq. ( |5q ) to compute this slope given 
a to almost equal to zero near the axis. 

The resulting system of equations is then of the form 



"'^ da 



^riij J — Ji 



,3 = 1-6 



(46) 



where A is its matrix. We use standard integrators for stiff 
systems of initial value problems of first order ordinary 
differential equations. We have chosen to study the prob- 
lem with initial conditions on the polar axis for simplicity. 
Other ways of numerically solving the problem could have 
been chosen. For example different boundary conditions 
could have been imposed for solutions that do not fill all 
space, arguing that thick disks would block some part of 
the available space. A pressure balance condition at the 
external boundary of the outflow could have been intro- 
duced. Another possibility would have been to prescribe 
densities or positions of the critical surfaces on the equa- 
tor. The solution on the equator would induce the solution 
on the polar axis. We have chosen to study systematically 
problems with initial conditions for numerical convenience 
but all the other possibilities have been tried and indeed 



converge to the same solution for similar parameters. This 
study of other boundary conditions for numerical integra- 
tion of the system has allowed us to check the numerical 
accuracy of the solutions. 

The first difficulty in solving this system consists in 
evaluating the initial values of the unknown functions on 
the axis. We return to this aspect later. We integrate 
this system for different values of the adjustable functions 
f2(a), Q{a) and of the parameter ao, the mass to magnetic 
flux ratio on the axis, which will be related a-posteriori to 
the total mass loss rate. The function Q{a), related to the 
entropy, is given by boundary conditions and depends on 
the object to be modeled. The rotation rate i^{a) of a mag- 
netic surface is very similar to that of the wind-emitting 
object at the base of the magnetic surface if this object 
is treated as a point source of wind. Once the positions 
and densities are computed at the three critical points it 
is possible to deduce all the other relevant variables such 
as the first integrals E, a, L or /3, w and the components 
of the magnetic field and of the velocity. 

Dimensional quantities The input parameters of the 
model can be selected so as to reproduce, at least qualita- 
tively, observed situations. Once given the mass Af* of the 
wind-enutting object, star or disk, the radius i?,, the tem- 
perature T,, the density rip, the total mass loss rate M,,, 
the magnetic field i?,, the factor Q* and 7, the dimension- 
less parameters O, Q, ao can be deduced. The parameter 
ao can be a posteriori related to the mass loss rate M*, 
i?*, and the magnetic field B*. So we define : 



{mpUpJ^ 
AttRIB^ 






m ■ 



(48) 



(49) 



All those quantities are nondimensionalized to reference 
values by setting Q = Q*/Qref, ao = at:/aref and fl = 
n^/riref- Major quantities of reference are given by 



Rref = 8.7 X 10"'m 
Pref = 3.4 X lO'^^kg.m' 
Pref = 5.2 X IQ-^Pa. 



(50) 
(51) 
(52) 



It has been found convenient to start by simply consider- 
ing constant values of Q{a) and Q{a). In the following sub- 
section the results of our model will be illustrated by con- 
sidering the specific examples of the Sun and of a normal 
TTauri star. We shall discuss in the next subsections the 
effects of varying rotation, mass loss rate and base tem- 
perature. Finally we shall show solutions for non-constant 
fl and Q. 
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3. The Solar wind (w < 1, fi == cste) 

Let us consider slow magnetic rotators, like the Solar wind, 
which we define as winds for which w <C 1. For slow ro- 
tators analytical solutions of the criticality equations can 
be found. We want to calculate the three first integrals 
which result from the imposition of regularity conditions. 
In particular, the specific energy is expressed as 



E{a) 




(53) 



When Ra is close to Rf we need to evaluate the ratios in 
this expression in a limit sense. Therefore to calculate E, 
we need to obtain the slope p at the Alfven point which 
is also a quantity of interest in the context of the Alfven 
regularity condition, Eq. (]l7|) . This slope (see Eq. (p^) can 
be expressed as 



p — 1 — l/zu. 
with 



A^ 



■t^tJ-OPA 



- iQpa' 



GM 
Ba 



2A2 



ft-Af^OPA 



n^cos^BR^, 



(54) 



(55) 



The specific energy can be evaluated only in terms of vari- 
ables defined at the Alfven point thanks to this variable 
■1x7. Then we get 



Eia) = 



A^ 



GM 



2^oPaR\ Ra 

+ ;; ^(1 + 



7 



7-1 



QpI 



(56) 



The specific energy is then fully determined in terms of 
the input parameters. 

Assuming a small value for the thermal parameter /3 
and eliminating e and g between equations (pi]),(p2), (RSh 
and (B4), we obtain equations for the positions and densi- 
ties of the slow and fast magneto-sonic critical points rel- 
atively to the positions and densities of the Alfven point. 
The position and density of the slow point exhibits inter- 
esting behaviors with respect to the different input pa- 
rameters. They are given in the cj <C 1 limit by 



rA 



PA 



/3(5-: 


57) 


2(l+/3- 


-C.2) 


2(l + /3- 


-C.2) 



4(^,-1) 



/3(7-l) 



1/4 



/3(5-37) 



(57) 



(58) 



In those equations the polytropic index 7 has two critical 
values, namely 7 = 1, the isothermal case, that does 
not create any problem and can be studied apart, and 



7 = 5/3. The latter value is never considered since 7 must 
be smaller than 3/2 for accelerated winds (Parker( [l963 ) ). 
When the thermal parameter /3 decreases, the slow point 
moves away from the Alfven point, approaching the source 
while the density at this point increases. The dependence 
on uj is weak, so the slow point is not much affected by the 
rotation. The fast point parameters present an opposite 
behavior. We have 



^ =1- 



4.^(1-^^ 



2 + w2 _ 2 



\2{\ 



/3(5-37) ^ "-^-^ fPi'f-i) 



+f3-uj^) 



PA 



4u;' 1 - 



= 1- 



/3(7-l) 



/3(5-37) 
2(l+/3-(.^2) 



O — df 
\ 4t-4 / 



/3(7-l) 



3/4 



3/4 



.(60) 



r(59) 



They depend much more on w than on (3. As lu in- 
creases, the fast point separates more and more from the 
Alfven point. When the rotation vanishes, the fast point 
approaches the Alfven point and both densities tend to 
become equal. Thanks to those results, analytical expres- 
sions can be derived for the dimensionless specific energy, 
e, and for the gravity parameter, g, namely 



l + (3 + 3uj^ 



9 



= 2^/4(7-1) 



3/4 



5-37 



2(l + /3-w2) 



/3 = 



(61) 



(62) 



While the relative specific energy e is almost unity and 
depends weakly on /3 and on w, the parameter g shows 
a clear relation with /3. Neglecting the effects of rotation 
and taking the cold limit, the specific energy becomes 



E ^ 



A^ 



2^oPaRa 



The gravitation parameter tends to vanish as 

5 — 3 7 
^3/4 I ii " ^ 



g^2^'^{^-lf 



37 



/32(-,-l) < 1. 



(63) 



(64) 



We need to compute the Alfven slope. We use Eq.(|6l|) in 
Eq.(p9|). We find a simple expression for the Alfven slope, 
valid in this low-rotation limit 



'^'-T2 



(65) 



All the previous results are expressed in terms of the po- 
sition and the density at the Alfven point. We now can 
calculate them. Using Eq. ( |65| ) and Eq. (|5|), we find a 
simple relation between the position and the density at 
the Alfven point 



PA 



[I + V2)A2 

iQR\ 



1/7 



(66) 
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The position of the Alfven point Ra can be obtained by 
integrating the Alfven regularity condition, Eq. dl^) which 
introduces a constant of integration Ci related to bound- 
ary conditions such as the total mass loss rate. We have 



Ra = {^j-qCipaE) 



2v/2 



(67) 



From Eq.(0), Eq.(p3h and Eq.(|67|) we can express Ra as 



Ra= ( 



CifioA' 



2(v^+2) 



(68) 



Thus in the cold limit it is possible to derive, from the 
initial conditions, all the relevant quantities, namely the 
positions and densities at the critical points, and the three 
first integrals of the motion, with the fast point almost 
coinciding with the Alfven point in the cold limit. So, the 
problem is fully solved. For slow rotators, we have found 
simple expressions for e, g, and for the Alfven radius and 
density with the introduction of a constant of integration 
given by initial conditions. 

3.1. Limit of vanishing rotation 

When the rotation vanishes the solution for Q{a) indepen- 
dent of a and a uniform flux distribution on the source 
must be spherically symmetric. Moreover the positions of 
the fast magnetosonic and Alfven points merge in the cold 
limit. We calculate the remaining four unknowns by mak- 
ing use of the slow point criticality equations. After some 
simple calculations, p^ and Rs are obtained as 



Ps = 



Rs 






(GM) 



7+1 



(2a2A2)7-i(2ag)2 



(69) 



(70) 



The Bernoulli equation written at the Alfven point gives 
an equation for Ra 



Ra{E-Q -p\ +^— = ^22 - 

7-1 Ra 2^0 a 



(71) 



where E can be found from the position and density at 
the slow point 



a'A^ GM 7 ^_i 

^RsPs Rs 7-1 

When P is small compared to unity, we get 
A' ^. Rtpl 



Ra 



2p,n^a^E n^a^' 



(72) 



(73) 



Using Eq. ( pSJ) , we can express in terms of the input param- 
eters the constant Ci which appeared in the integration 
of the Alfven regularity relation 



Ci = 



2 /i?4^2x2(^+^) 



HoA'^ VMo"^ 



(74) 



Using the expressions of Rs and Ps given by Eq. ( |70| ) and 
Eq.(p9[), this can be also expressed as 



Ci 



^2^9+4V2^8(v^+2) y 2''-lQ5a27-l 



„ 8(%/2 + 2) 

ji{GM)''-^A^-^''\^^^ 



(75) 



In the case of vanishing rotation, it has thus been 
possible to obtain an analytical solution. We got simple 
expressions for the positions and densities at the critical 
points and the constant of integration in terms of the in- 
put parameters only, this including the wind mass loss 
rate, represented by parameter ag- This solution is im- 
portant because it gives initial conditions on the axis for 
the numerical integration of the equations of our model of 
rotating MHD winds. 

3.2. Slow cold rotator 

Ignoring /3 in Eq. (^ and Eq. (|6^) we calculate the solu- 
tion for the fast point parameters 



Pf = Pa\1 
r/ = r^ ( 1 - 



4cj2 



(2 + ^2)2 

_JuP_ 

{2 + LU^f 



Pa{1-lo^) 



TA 1 



(76) 



(77) 



For w <C 1, the relation gives rise to finite values of the 
quantities. The fast point remains close to the Alfven point 
and depends on uP' . 

3.3. A numerical solution for the Solar wind 

Taking into account realistic values for the sun it has been 
found that lo was equal to 0.55 with Q = 2.1, fi = 3.1 and 
a = 1.7 as input parameters. Thus the Sun can be consid- 
ered as a slow rotator with respect to our classification, or 
at least at the border between slow and intermediate rota- 
tors. The numerical solution shows clearly that the critical 
surfaces do not differ much from spherical shape, partic- 
ularly for the slow surface. The fast magnetosonic point 
is close to the Alfven point even near the equator, where 
the effect of the rotation is more important as expected. 
We have also plotted the various densities in Fig. g. The 
fluid near the polar axis is denser than near the equatorial 
part and one should keep in mind the large values of the 
densities for later comparison with fast rotators. 



4. Jets from YSOs ( (|j 



4.1. Stellar winds (H — cste) 



a; <C 1 



Now, we focus on the case of fast rigid rotators. Fast ro- 
tators are defined here as having a rotation parameter w 

close to (I) ^ ■ This value is a maximum for uj in this con- 
ical model and must be regarded as an asymptotic value 



10 



T.Lery et al.; Outflows from magnetic rotators. I 



U.J 


















^^^~~~— ~„_;;~ ^ 






0.4 


- 






^\. ""^ ^ 


- 


S 0-3 


- 








N 


W) 










\ X 


OJ 










\ \ 


■^ 0.2 


— - 


■\ 


\ 




\ \ 
\ ^ 


0.1 
n 






\ 
\ 
\ 

1 1 


1 


\ \ 
\ \ 
\ \ 

I \ 

II 1 1 



0.2 0.4 0.6 

equatorial radius 



0.8 




Fig. 2. results obtained for defined by Q = 2.1, 51 = 3.1 
and a = 1.7 as input parameters. Quantities referring to 
tiie Alfven point are plotted with solid line while dashed 
lines correspond to the fast point and long dashed line to 
the slow point. 

for plasmas rotating at a very large physical rotation rate 
n. Similarly to the case of the slow rotator we can de- 
rive analytically in this limit the positions and densities 
at the critical points relative to those at the Alfven point. 
Assuming that r^ <C r^ and ps ^ pA they are given by 



TA 


\ 2 ) 


//5-37 


Ps 
PA 


V 2 V 


--"))" 



3w3 -2w 



A \ 4(t^1) 



(79) 



The slow point parameters depend both on w and (3. 
When /3 decreases, if the Alfvenic point remains approx- 
imately at the same position and density, the slow point 
approaches the source and the corresponding density in- 
creases. At the fast point assuming that rf ;:^ ta and 
Pf ^ Pa we have similarly 



rA 



1 + 



4 
3a; ^ 



{^^y 



'^(2-l)V<^-" 



EL 

PA 



'5-37 

7-1 



(rA \2 



(iuj^ - 2lj'^\ 



( 7 + 1 'l 

V4(t-1); 






2/3N 



(80) 



(81) 



If (3 approaches zero, the fast point is rejected to infinity, 
a well known result in the cold plasma limit (Kennel et 
al. (1989)) The dimensionless specific energy e can be de- 
rived by substituting these expressions for the positions 
and densities at the critical points in the Bernoulli equa- 
tion. We obtain: 



3^3 



TA 



Lj^n 



and if rf ^ va it can be reduced to 

4 

e = 3aj 3 . 



■ 6a; 



2/3 



(82) 



(83) 



This gives the specific energy i? as a function of the density 
and position of the Alfven point as: 



E 



2V^ 



(84) 



Using the definitions of uj and g which in the fast rotator 
limit can be written a,s g — 2uP' [r s / r Ay^ , the position of 
the slow and the fast points are found to be given by 



GM 



f72cos6l2 



1/3 



(85) 



As the rotation increases, the slow point gets closer to the 
source in proportion to 51^^/3 ^j^g fg^g^ point position is 



2^4 



a^r 



rf 



GM 



(86) 



since ta will be found to increase with the rotation rate, 
this expression shows that since va will be found to in- 
crease with the rotation rate, the fast point is rejected far 
from the Alfven point when the rotation grows very large. 
Interesting physical consequences will be given in the final 
'(TOcussion part. 

4.1.1. Fast cold rotator 

If moreover the wind is cold, the fast point goes to infinity, 
in the limit of vanishing entropy. 



(87) 



lim Of = 
lim rr — oo 

/3^0 ■' 



It could be feared that the geometry of magnetic sur- 
faces upstream from the fast point might not be coni- 
cal at very high rotation rates. Numerical exploration of 
the positions of critical points in models with non-conical 
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geometry having a cylindrical asymptotic shape has how- 
ever shown that, provided the rotation does not grow to 
extreme values, the fast point still remains in the quasi- 
conical region surrounding the source and does not shift to 
the asymptotically cylindrical region. Therefore our sim- 
ple assumption may have a somewhat broader domain of 
validity than would be anticipated. This is why we felt 
that it is still worth pursuing the study of the approxi- 
mate model even in this limit where its validity becomes 
disputable. When (3 and g are small compared to unity, 
we find 



w 



(89) 



The Alfven slope, for the fast rotator regime, approaches 
a constant value 



(90) 



We can also deduce the radius and the velocity at the 
Alfven point as well as the specific energy and angular 
momentum 



Ra^ 



VPA 



1/2 



A 



1/3 



2) \iJLQaQ.cos(j 
A 2 fAn'^cos^e'^^''^ 



/ioai? 
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2 \ /ioct 
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L 



3 f An^/^cos'^e 
2 V ^0" 



(91) 
(92) 
(93) 
(94) 



we thus obtain expressions of the integrals of the motion 
E and L relevant to this case. Inserting these results in 
the Alfven regularity condition and integrating it, we find 



16 /IT 

D 27 V 19 
Ha 



s6 



Ra 



2(1- 



Ea 



= C2 



(95) 



where C2 is an integration constant. This gives the Alfven 
radius in the form 



Ra 



Con 



Ea{cos( 



^2(1- 



'm 



24-2V5T3 



(96) 



Comparing Eq.(^l|) and Eq.(pl), we find that the constant 
C2 is related to other quantities by 



C, ^an 



3/2 



A 



2) noaQcosO 
When 9 — Q, it becomes 
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(97) 



(98) 



and 

Ea 

~n 



3 fanA^ 



1/3 



(99) 



This last relation is directly related to the escaping 
poloidal electric current. So we have obtained a simple 
expression for the net current. We have postulated that 
/? <C 1 in this part. This is express by the following in- 
equality : 



a< 



7-1 {Afl^ cosO^) 



2^,2/3 



37-1/3 
Mo 



(100) 



It shows that in the limit of our assumptions a is restricted 
for given value of all the other input parameters. We will 
come to this result later in the section dealing with nu- 
merical results. ^ 

For very fast rotators, i.e. oj ~ (|) ^ , we can easily 
express this constant of integration, which then allows to 
calculate all the variables and integrals of motion of the 
system. In this case the Eq. (^ closes the set of equations 
that define the solutions, as did Eq.(|75|) in the case of no 
rotation. 



4.1.2. An example of TTauri star: BP Tau 

This section is designed to show by an example the be- 
havior of the six variables for which our model gives solu- 
tions in terms of boundary conditions and of the deduced 
quantities. As typical, we present the results for a typical 

TTauri star, BP Tau. 

From Bertout et al. (1988), we take the reference val- 
ues as M, 



2 X 10^'^M, 



'-iQyr 



M* = O.8M0, i?, = 3R, 



■0: 



T* = 9 X 10^, Up = lO^CTO^^^ and B* ^lOOOG. We_deduce 
the dimensionless input parameters Q — 0.05, Q = 1.8 
and cJo = 0.1. The three critical surfaces (upper panel) 
and corresponding densities (lower panel) are represented 
on Fig. g. The lower part of the figure shows clearly the 
decreasing trend of the densities from the axis to the equa- 
tor and stresses the differences in magnitude of the den- 
sities at the different critical points, the density at the 
slow point being one hundred times larger than the oth- 
ers. This refiects the position of the critical surfaces. On 
the upper part of figure 2, it can be seen that the fast and 
Alfven critical surfaces definitely exhibit shapes different 
from spherical, particularly the fast one. One can notice 
that the further we go away from the axis the more the 
fast and the Alfven surfaces get apart. As seen in Sakurai 



(1987) and Belcher & McGregor (1976), as one gets closer 
to the equator the critical surfaces get more and more 
elongated along the equator. However a difference appears 
near the axis of rotation. In the latter cases, the critical 
surfaces get elongated along the pole. This is mainly due 
to the constancy of the angular velocity r2(a) and of the 
entropy Q{a). When this assumption is relaxed, equivalent 
behaviors have be obtained. On the other hand rotation 
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radius 




Fig. 3. Example of the results obtained for a typical 
TTauri star defined by Q = 0.05, Cl = 1.8 and a = 0.1 
as input parameters (7 = 1.2). In the upper graph, the 
densities at the three critical points are represented with 
respect to the magnetic flux a in linear-logarithmic frame. 
The lower graph shows the corresponding critical surfaces 
in dimension- less quantities. The Alfven surface is the 
solid line curve and the fast point surface the dashed line 
one, while the slow surface has a heavy solid line and is 
situated close to the origin. Note that, on the whole, the 
slow surface is almost circular while others are not at all. 



70" r 
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Fig. 4. Example of variations of the three first integrals of 
the motio n for a typical TTauri star (BP TAU, Bertout et 
al. ( 1988 )) with respect to the magnetic flux a. they are 
deduced from the positions and densities of the critical 
points. The solid line represents the specific energy E, the 
dashed line the angular momentum L and the long-dashed 
line the mass loss to magnetic flux ratio a. 



has little effect on the slow mode critical surface which al- 
most keeps a spherical shape. In fact the slow surface gets 
closer as the velocity increases, as one would expect for 
such a type of wave considering the increase of centrifugal 
acceleration near the equator. The behavior of this surface 
will be studied in more details in next sections. From the 
positions and densities of the critical points we deduce the 
other interesting variables of the problem, such as the first 
integrals of the motion E, L and a as represented on Fig. 
^. The specific energy seems to vary little as compared 
to the other quantities. Its variation is less than one per- 
cent of the mean value. Thereafter the amount of energy 
is slowly increasing from the value on the axis, given by 
Eq.(Q), to an higher value on the equator. The specific 
angular momentum increases from the axis and then re- 
duces its growth to level off at a value close to unity. The 
behavior of a is identical to that of the Alfvenic density 
as seen in Eq. (g) . One should keep in mind that the value 
of a on the axis is prescribed as an input parameter. 

In this subsection, we have shown a particular set of 
numerical solutions of the equations of this model. In the 
following sections we comment on the influences of the 
input parameters on the solutions. 

4.1.3. Effects of the amplitude of the Mass Loss 



2 - 



1 - 



-1 - 



-2 - 




Fig. 5. The Alfven surfaces for increasing ao are shown 
with solid lines. Circular shapes are represented with 
dashed lines for comparisons with the solutions. Smaller 
surfaces correspond to higher values of ao while most dis- 
torted surfaces have smaller ag- The smaller ao, the bigger 
the distortion and the higher the height of the Alfven sur- 
face on the polar axis. Small ao can either correspond to 
low mass loss rate, or to high magnetic flux. 
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Fig. 6. As shown on the left panel, /3 is small compared to 
one for various values of ao- The upper long-dashed curve 
corresponds to high ao while going down the curves, one 
finds solutions for smaller ratios. Those curves strengthen 
the assumption for our analytical calculations where /3 is 
considered small compared to one. On the right panel, the 
curve indicates the total mass loss rate against ao ■ Notice 
that the slope of the curve is almost unity. 



The mass loss rate M is a quantity that can be deduced 
from observations. In our model, it is represented by the 
parameter ao, the mass flux to magnetic flux ratio on the 
polar axis. The relation between ao and M can be worked 
out a posteriori. We have calculated solutions for different 
ao's, keeping for the other parameters the same value as 
in the previous calculations. The corresponding Alfvenic 
surfaces are presented in Fig. |[ Smaller surfaces corre- 
spond to high ao's while the distorted ones correspond to 
smaller values. This variation is similar for all the critical 
surfaces. The critical surfaces strongly inflate for decreas- 
ing ao and lose the initially spherical shape that they have 
at high mass loss rates. It shows that for strongly mag- 
netized systems, the structure of the zone of acceleration 
gets complex. 

Fig. M shows the evolution of /3 as a function of the 
magnetic flux for different ao's. Fig. p^ shows the simple 
relation which exists between ao and the mass loss rate 
with the other input parameters fixed. In Fig. o the values 
of (3, given by Eq.(p^), are represented against the rela- 
tive magnetic flux for a series of decreasing values of ao . 
When ao is large, the behavior of /3 is nearly linear but 
for lower values the curves start to decrease steeply and 
reach a low equatorial value. This strengthens one of the 
assumptions made to derive our analytical results, namely 
that the value of /3 is small compared to unity. The value 
of /3 on the axis depends directly on Q. If the entropy 
is large enough, (3 can become large too and even reach 



unity. This would be the case for very hot outflows. The 
smaller ao , the better the assumption of smallness of P 
is justified. Note that the critical surfaces for magnetized 
rotators with small mass loss rates are strongly inflated 
equator-wards and are in no way similar to the quasi- 
spherical shape they have for rotators with large M. In 
this section we have shown that fast and slow rotators 
developed different shapes of the critical surfaces. It does 
not seem possible to infer the solution for the fast rotator 
from the one obtained for a slow rotator by simple scaling 
scaling arguments as attempted by Shu et al. ( 1994 ). 



- 4.1.4. Rotational effects 




Fig. 7. Effects of the angular velocity fl on fast sur- 
faces with parameters of a typical TTauri star. Circles 
are shown to visualize the distortions of the critical sur- 
faces. These curves indicates that the faster the rotation, 
the bigger the distortion with respect to sphericity is. The 
position of the critical points on the polar axis is of course 
independent of the rotation. 



The value of the angular velocity has a major influence 
on the behaviors of the solution. In Fig. ^ several Alfvenic 
surfaces are represented for different il. One should notice 
that the central emitting object is reduced to a point and 
does not coincide with the first circular surface. This inner- 
most surface corresponds to the smallest rotation rate and 
the other surfaces expand off the central part for larger 
and larger rotation rates. The faster the rotator, the larger 
the change in the shape of alfvenic surfaces. The effect of 
the rotation is strongest at the equator. The Alfven sur- 
face remains at a finite distance from the source for finite 
rotation rates. The Alfven point only goes asymptotically 
to infinity for ultra-fast rotators. 
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Fig. 8. The variation of the rotation parameter to with 
to the magnetic flux a is plotted for different values of the 
rotation rate fi . The asymptotic limit of (|)2 is repre- 
sented with an heavy dashed line. Lower curves represent 
slow rotators while upper curves stand for fast rotators. 
An arbitrary limit separate the two type of rotators at 
w ~ 1. Close to the axis, a slow rotator type behavior 
always remain whatever the rotation rate. 



slow rotator 
Fast rotator 




In Fig. g we present a plot of Ci;(a) for different rotation 
speeds with the same previous input parameters Q and ao- 
The lower curves illustrates slow rotation rates and those 
above them correspond to increasingly larger rotation. 
The horizontal dashed lines represents the asymptotic 



limit for ultra-fast rotators where co = 3/2' 



3/2 



1.837. 




The curves pass smoothly from the region of slow to fast 
rotation. Once the asymptotic value for uj is nearly reached 
on the equator, the asymptotic region extends towards the 
axis for larger rotation reducing the extent of the region 
where the rotation parameter stays small. This region re- 
gion never disappears completely, though, because obvi- 
ously the polar axis itself is necessarily in the slow rotation 
regime. 

For fast rotators the Alfven point is at a large dis- 
tance. In this case, we see that the alfven speed is large 
compared to the sound speed. At the equator the slow 
magneto-sonic surface as seen in Fig. g gets closer to the 
source as the rotation rate increases. The effect of the cen- 
trifugal acceleration is more and more pronounced as one 
reaches the equator. The slow mode speed gets closer to 
the sound speed, and the slow mode acquires the charac- 
ter of a sound wave guided along the field line. For this 
reason we can see the slow point on the equator moving 
to the source. The radius of the slow mode surface on 
the equator approaches the radius where the sound speed 
corotates with the Keplerian speed, given by: 



fGM\ 



1/3 



(101) 



We have represented the radius at the slow point situated 
on the equatorial plane as a function of increasing rotation 
speeds on the right panel of Fig. g together with the pure 
sonic point. It is clear that the slow magneto-sonic point 
tends to the sonic point. We now have a view of how all 
the critical surfaces behave as the rotation rate vary. Fast 
magneto-sonic and Alfven surfaces strongly inflate when 
rotation increases while the slow surface gets narrower to 
converge to a nearly cylindrical shape. It seems that an 
increase of ft has almost the same effects as a decrease of 
a, except that the positions of the critical points do not 
change on the polar axis when the rotation rate varies. 



equatorial radius 

Fig. 9. The shapes of the slow magneto-sonic surfaces for 
increasing rotation rates Vt are shown on the left panel. 
The slow surface has a opposite behavior compared to the 
fast and the Alfven surfaces. It deflates while il increases. 
On the right panel of the figure, the numerical solution for 
the magnetosonic slow point is represented with the vari- 
ations of the pure sonic point with respect to the rotation 
rate. This shows that the slow magnetosonic point tends 
to the trend of the sonic point for increasing rotations. 
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4.1.5. Thermal effects 

We have assumed a polytropic equation of state. Thus the 
two parameters that we can vary to illustrate the thermal 
effects are the entropy and the polytropic index. We con- 
sider these parameters as constants as a function of flux 
variable a. In these calculations, $7 and ao are fixed. When 
Q, the factor related to the entropy, decreases, the critical 
surfaces inflate, the rotation parameter uo increases and a 
reduces. This can be understood as follows: the more the 
plasma is heated, the lesser is the influence of the rotation. 
When the polytropic index is varied and the thermody- 
namics turns progressively from isothermal to adiabatic. 
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Fig. 10. Some of the thermal effects are represented on 
this figure. The upper panels show the variations with Q 
and 7 of the distance to the origin of the critical points in 
the equatorial plan. Below we show how to and a vary, 
smaller entropies make the rotation parameter smaller 
while they produce larger ratios of mass loss rate to mag- 
netic flux (Note that the entropy decreases from left to 
right). An increase of 7 have the same properties. The 
isothermal case produces the larger mass loss rate and as 
7 increases the rotator turns to be faster. 



the influence of rotation becomes more and more impor- 
tant. This is because the effect of gas pressure on the dy- 
namics is maximized by an infinite thermal conduction. 
For most solutions the ratio of the fast mode radius to 
the Alfven radius is not very large. Thus the variations of 
Q and 7 have a non-negligible influence but nevertheless 
weaker than that of varying Q and ao . A reduction of the 
entropy has a similar effect to that of an increase of the 
polytropic index 7. The smaller the entropy, the closer the 
flow becomes to that of a fast rotator. 



4.1.6. Limit cases 

In the previous section, we have seen that the fast ro- 
tator regime is achieved before Q vanishes. Similarly, all 
the other parameters vary in a restricted range. For ex- 
ample, as the rotation Q increases, the rotation parameter 
to reaches the maximum value (l)^ . Numerically we find 
that there is an lower limit to ao for given il and other 
parameters. To understand it, let us recall the equation 
g = 2GMfila^R\A-'^. 

It has been found that 



VlcosO iiQaR\ 



So for fast rotators the Alfvcnic radius tends to be almost 
equal to 



1/2 



( A \ 



1/3 



The gravity parameter the becomes 



2GM 

TAVPA 



3/2 



l'2GMfio 
V Acos9'^ 



a 

n 



(103) 



(104) 



Meanwhile we have found that this parameter could be 
expressed as 



9 



s 4 * 



(105) 



It shows that g is small with respect to unity and that 
it does not vary much since Xs mainly depends on the 
thermal parameters and not on much on the mass loss 
rate or the rotation rate in the case of fast rotators. It 
also shows that in the case of very fast rotators either the 
angular velocity must be very important or the mass loss 
rate must be small to insure that g is small. One can also 
deduce that there is a maximum Q for a given such as a 



nr^ 



^maxf{(Xa,Q) 



(106) 



The fitting of the numerical limits of our model where 
u! — (I) ^ agrees with this result and gives 



nr. 



1 



(107) 



(102) 



Thus this shows that if the outflow is to support a given 
mass loss rate, it cannot rotate at any large rate. The 
converse is also true: if the magnetic field lines turn at a 
given rotation rate O, there is a minimum mass loss rate 
from the central object. 

4.1.7. Discussion 

We have presented an extended set of solutions of the 
model, for constant Q and Q across the magnetic field 
lines. We obtain the first integrals of the motion, the com- 
ponents of the velocity, of the magnetic field and all the 
variables characterizing the flow. This study revealed some 
interesting properties. We have shown that an increase of 
fl up to the fast rotator regime has a similar effect to a 
decrease of ao, or Q, and to an increase of 7. Fast rota- 
tors display critical surfaces which are strongly inflated 
equator- wards, that drastically differ from their quasi- 
spherical shape in the slow rotator regime. The faster the 
rotator, the more valid the assumption that P is small. 
The influence of thermal parameters is smaller than that 
of rotation or mass loss rate but not negligible. Finally we 
have also found that a given rotation rate rate compels 
the mass loss rate of the outflow to have an inferior limit. 
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Fig. 11. The variations of the rotation rate Q. with a for 
the differential rotators, considered in this work. The dot- 
dashed Hne represents a Keplerian rotation. After from 
the bottom, one find a step-shape rotation with Vt varying 
from 1 to 0, a step-shape rotation from 1 to 0.5, a decreas- 
ing stellar type rotation, a constant one, a normal stellar 
type differential rotation and finally a step-shape rotation 
from 1 to 1.5. 



4-. 2. Other variations of fl and Q 

The distribution of outflowing material in space and ve- 
locity is not well determined observationally, but some 
properties of such flows have still be recognized. The rel- 
atively wide spectral lines which have been measured at 
most places in such outflows indicate that material with a 
range of velocities is present on the line of sight. However, 
in a number of well-collimated outflows, such as NGC2024 
and NGC2264G, there appears to be a shear flow, with 
higher velocities near the polar axis, and lower velocities 
at t he sid es of the flow (Margulis et al. ( 1990 ); Richer et 
al. ( |1992| )). These observations motivated us to look for 
outflows with rotation rates and entropy varying with the 
flux variable a, that could reproduce such behaviors. 



4.2.1. Non-constant Q 

We now describe some specific forms of the rotation rates 
and discuss trends and properties that can be identified. 
For the rigid body rotation, the angular velocity is con- 
stant with flux a. The stellar-type rotation is inspired 
by the equation of differential rotation of the sun. We 
have modeled a stellar-type rotation with angular veloci- 
ties larger at the equator than at the pole. We have also 
formulated variations of the rotation rate that could ac- 
count for t wo co mponents outflows such as the models of 
Shu et al (1994). A central fast flow is surrounded by a 
slower wind. We present two different differential rotators 



Fig. 12. Fast critical surfaces for a differential rotator. 
In this example the rotation rate starts with a fast rota- 
tion on the axis and then reduces to a smaller rotation 
rate on the equator. The surfaces start just like fast rota- 
tors shown on Fig.(0) but change drastically their shapes 
afterwards giving this butterfly shape. 

with this characteristics. We call them two- velocities out- 
flows with a step-shape rotation. The analytical expres- 
sions of all these differential rotation rates can be found 
in appendix C. We have plotted them as functions of the 
magnetic flux on Fig. 11, Solutions of the problem with 
such profiles of rotation differ from previous ones. Criti- 
cal surfaces, densities, and all the other variables, strongly 
vary with a. Fig. |lj is an illustration of of the fast crit- 
ical surfaces for a two-velocities outflow. This produces 
"butterfly" shapes showing the possibility for two kinds 
of outflow with different characteristic to live side by side. 
We have also plotted the solution for a central rigid 
rotator surrounded by a keplerian disk in Fig. O. The 
central part near the axis of rotation has got a constant 
rotation velocity up to a = 0.3, then fl follows a keple- 
rian variation. This kind of model could apply to Young 
Stellars Objects where the central object is launching an 
0-win d and the accretion disk an X-wind (see Shu et al. 



( 119941) ) 



4.2.2. Non-constant Q 

The variation of the entropy with magnetic flux a seems to 
have little influence on the solution in our computations, 
presumably because /3 remained small either due to small 
Q values or to the rapid rotator regime we imposed in 
several cases. Nevertheless we have studied different pos- 
sibilities and only in in limit cases of the value of Q and 7 
the results change. In previous sections we have presented 
solutions for constant entropy with different values. We 
have also tested different variations with a of the func- 
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At the fast point the solution for these two parameters is 
given by 



6 = 3C.4/3 
and 

2 
.9= — 



^^ + 3^2/3 _ -4/3 



(110) 



(111) 



It is possible to deduce equations for the positions and 
densities at the slow point if one considers that j/s 3> yA 
and Xs <C XA and gets 



" ^2(1+^23.2) 



(112) 



TA 



//?(7- 1)^(^-1) 



5-37 



(7-l)(3a;4/3_2w2) 



(113) 



Using the last two equations one gets a better equation 
for the density 
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Fig. 13. Positions and densities of the the slow, Alfvenic 
and fast magnetosonic critical points for a fast rigid rota- 
tor surrounded by a keplerian disk. 



tion (5(a), with increasing or decreasing entropies from 
pole to equator and a step-shape profile. Details can be 
found in appendix C. The general trend is in line with 
the results obtained in our study of thermal effects. For 
either increasing or decreasing entropies, critical surfaces 
respectively get closer to the source or inflate outwards. 
For step shape variations the critical surfaces are similar 
to surfaces represented in Fig. ^. 

5. Intermediate rotator 

5.1. Analytical results 

Considering lu close to unity and neglecting /3 with respect 
to unity it is again possible to obtain analytical results for 
the energy and gravitation parameters. At the slow point 
we get 



■Py. 



,7-1 



4 ? 



2lu^ 



and 



2^3 



g = 2uj^x' 



(108) 



(109) 



£s_ 
PA 



/3(5 - 37) 
2(3cj4/3 _ 2^2) 



(114) 



It is also possible to get the variations of the position 
and density at the fast critical point. In the intermediate 
rotator regime, no direct simplification can be done on the 
positions of the critical surface, so the full set of solutions 
is hard to find analytically but nevertheless the solution 
is fully determined. Only numerical studies can help to 
understand the behavior of such solutions in this regime. 
This is what the next part intends to do with comparisons 
of the numerical and analytical results. 

5.2. Analytical vs numerical results 

The analytical solutions and the numerical results ob- 
tained for the specific energy and presented in the pre- 
vious parts have been plotted in Fig.(|lJ) for comparison. 
The figure shows the specific energy against the rotation 
parameter lu. The results agree in the slow regime and in 
the very fast regime, as expected. The intermediate zone 
that we have arbitrarily situated between lu = 0.75 and 
LU = 1.5 shows clearly the limit of validity of the solutions 
for the slow and fast rotators. It shows a smooth transi- 
tion between the different categories of rotators. In fact, 
this region is narrow in the physical space for a fast ro- 
tator since it is the transition between the central slow 
part of the jet close the axis of rotation and its fast part. 
So the analytical results are a good approximation of the 
numerical solution. Similar comparisons have been made 
for the positions and densities of the critical points with 
similarly positive conclusions. 
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Fig. 14. Comparisons of numerical results and analyt- 
ical calculations of the specific energy. On the left, e is 
represented as a function of the magnetic rotator energy 
u>. There are good agreements for slow and very fast rota- 
tors. This is not the case of the intermediate region around 
w = 1 that has no analytical solution. But this region is 
really narrow with respect to the magnetic flux a as shown 
on the plot on the right hand side. 



6. Discussions 

6.1. Comparison with other models 

In a review of the theory of magnetically accelerated out- 
flows and jets from accretion disks, Spruit ( |1994D discusses 
how the wind structure should depend on mass flux (see 
also Cao fc Sp ruit ( 1994 ) ). Using the model of Weber 
& Davis ( |1967D he concludes that when the mass flux is 
small, the Alfven radius extends far away from the ori- 
gin. Our analytical and numerical results agree with these 
conclusions, with even more generality, since our model 
fllls all space and does not assume any a-priori variations 
of the different variables. We find, using the conditions of 
vanishing of the differential form of the Bernoulli equation, 
that, if /3 is small compared to unity, the Alfven radius on 
the equator can be approximated by 



Ra 



-1/3 



(115) 



For large mass loss rate. Spruit finds that the Alfven point 
does not recede arbitrarily close to the origin but reaches 
a minimum value. In our model when a becomes large 
the fast point can be as close as possible to t he so urce. 
Following the formalism described by Spruit ( 1994 ) and 




used by Sakurai (1985) we can deduce the relation between 
LJ and a that is relevant to fast rotators, which is 



(116) 



This solution is represented in Fig.(|15D together with our 
numerical solutions on the equator. When a decreases, 
the two curves get closer while rotation parameter grows 
larger and reaches (3/2)^/^ when a = 0. The ultra-fast 
rotator limit is clear in this case. Again we find that the 
larger the mass loss rate the smaller the rotation parame- 
ter. A disagreement arises with the previous formula, since 
it does not allow a to be larger than unity. This shows the 
limited region of validity of this analytical approximation, 
which nevertheless describes quite well the very fast rota- 
tor. Our analytical solutions tend to Sakurai's results in 
the cold limit. In our case the slow point is defined by 
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1/3 



This is equivalent to 
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and the specific energy is 
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Moreover, we can express those quantities as functions of 
our input variables for very fast rotators and we get 



GM /Acose^^^^"^ 
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and 
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(120) 
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Thus we find the same results as those described by Spruit 
in the cold limit, with more generality since we solve the 
transversal force balance equation on the Alfven surface. 
The shapes that we have found for the critical surfaces 
differ from those obtained or assumed in many other mod- 



els. For example, Sakurai (1987) found ellipsoidal shapes 
with his numerical studies. Blandford & Payne ( |1982 ) have 
chosen conical Alfvenic surfaces in their self-similar model. 



Chan and Henriksen (198C) and Sauty (1993) used flat sur- 
faces perpendicular to the axis of rotation. Our solutions 
show rather different shapes. It appears to be difficult to 
capture realistic behaviors by modeling the critical sur- 
faces with simple functions. 

Further developments within the framework of our 
model could be done to investigate the evolution of the 
angular momentum of solar-mass stars during pre-main 
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Fig. 15. Comparisons of numerical results and analyti- 
cal calculations of the relation between uj and a on the 
equator. The upper curve (solid line) represents the nu- 
merical result while the lower one (dashed line) the ana- 
lytical solution as derive following the formalism of Spruit 
( |l994| ). The entropy is taken as constant. The two curves 
converge to each other when the rotation parameters be- 
comes larger, which happens when a decreases. We find 
again that the bigger the mass loss rate the slower the 
rotation parameter. 
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Fig. 16. Comparisons of numerical results and analyti- 
cal calculations of the relation between the Alfven radius 
on the equator and lo. The upper curve (solid line) rep- 
resents the numerical result while the lower one (dashed 
line) the analytical solution. The two curves are similar al- 
most everywhere and particularly for very slow and very 
fast rotators. 



and main sequence phases. Charbonneau (1992) recently 
derived a braking rate from a Weber & Davis (1967) model 
of the solar-wind assuming a dynamo relationship where 
the stellar magnetic field scales linearly with the stel- 
lar angular velocity Jl*. He finds that the braking rate, 
dfl/dt, asymptotically scales as fll at low fJ*, and as Ql 
at high n^,. In the Weber & Davis model, the flatter Jl*- 
dependency of the braking rate for rapid rotators results 
fro m the change of the magnetic wind structure (Mes- 
tel (|196S|) , Belcher & MacGregor (|l976|) ). That distinct 
braking laws apply to slow and fast rotators is further 
supported by observations: Skumanich's ( 1972 ) relation- 
ship (ri* ~ i^^/^), which follows from an fi^ braking law, 
is valid for slowly rotating dwarfs but fails for younger, 
more rapidly rotating zero-age main sequence star whose 
angular evolution is more appropriately described by an 
nl braking law (Mermilliod & Mayor ( |l990[ )). Thus our 
work takes place in a context of interesting studies and 
some further investigations should be done in this direc- 
tion. 



6.2. A criterion for the classification of rotators 

With respect to the previous results, it has been interest- 
ing to neglect the effects of the thermal parameters and 
only to look at those of the rotation, of the magnetic field 
and of the mass loss rate. It has been shown that an in- 
crease of ft has effects similar to those of a decrease of ao , 
all the other parameters being fixed. This suggests that 
the ratio — might be relevant for a study of magnetized 
rotators. That is the reason why we have plotted the rela- 
tion between the ratio fi/ao and oj. First this gives a global 
relation between the input parameters of the model that 
can be related to physical parameters given by observa- 
tions (See appendix C) and the variable to that is a result 
of the calculation. 

We find a continuity between slow and fast rotators. 
Slow rotators correspond to w <C 1. Since this class con- 
tains the Solar wind, winds of of B[e] stars, AGE stars, it 
can be identified with winds in general. In the intermedi- 
ate regime that arbitrarily corresponds to 0.75 < tu < 1.5 
we find both O stars and TTauri stars. This region clearly 
shows a smooth transition between slow and fast rota- 
tors. In the fast rotator zone where uj is close to (|) 
only object with jets, i.e. outflows of TTauri stars, of cata- 
clysmic variable, of microquasars such as SS433 and of ac- 
tive galactic nuclei, are to be found. This shows two main 
types of objects, winds and jets with a transition in be- 
tween. One should notice an important difference between 
the TTauri stars belonging to the intermediate zone and 
those in the fast rotator regime. They have larger mass 
loss rate (> lO~^M0yr~^) compared to the other ones 
(ilO~^M0yr^^). This could be due to a difference of ages 
for these TTauri stars. This suggests that M and oj are 
closely related. Then our model could provide a tool to see 
the evolution of this type of object and to classify them. 
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Fig. 17. Classification of the magnetic rotators, the log- 
arithm of the ratio — given by the observational lit- 
erature is plotted against the rotation parameter uj ob- 
tained by the calculation for several astrophysical objects. 
The grey zone represents the region of uncertainty over 
the input parameters. The two vertical thin dashed line 
shows the approximate limits between slow, intermediate 
and fast rotators. The limit for ultra-fast rotators is at 
(1)3/2 « 1.837. 
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From the input parameters related to astrophysical 
quantities given by observations and with the results of 
the model, we can classify outflows of magnetized rota- 
tors into slow, intermediate and fast rotators. We have 
also computed solutions for planetary nebulae and Wolf- 
Rayet stars but, since the origin of their wind is mainly 
radiative, they should not be compared with the other 
objects. Nevertheless the numerical results obtained for 
these types of astronomical objects shows that they are 
definitely slow rotators following our classification. 

6. 3. Limits of the model 

Our model possesses several advantages. First it is not 
self-similar and secondly it can accept any boundary con- 
ditions. Moreover the simplifications allows us to find fully 
analytical solutions in limit cases and systems of equations 
which are easy to solve numerically in all cases. The model 
is not simply a continuous set of closely packed Weber- 
Davies solutions but it is made coherent with the Alfven 
regularity condition on the Alfven surface, which retains 
part of the cross-field balance. It also gives the possibility 
to solve the Bernoulli equation from the source to the fast 
point and finally the asymptotic resolution of the outflow 
is possible (See Lery et al. ( |1998b[ )) giving the possibility 
to relate the properties of the source to the asymptotic 
behavior of the outflow. 



The main weaknesses however are the following. First 
our model does not give an exact solution to the system 
of coupled Bernoulli and transfield equations. Secondly 
the geometry of the magnetic field lines is idealized. This 
latter difficulty could be avoided by iterating on the field 
geometry. 

7. Conclusions 

In this Paper, we have proposed a simplified set of the 
equations for rotating magnetized outflows by assuming 
the shape of the poloidal magnetic field lines up to the 
fast magneto-sonic point. Rather than solving the equilib- 
rium perpendicular to the flux surfaces everywhere, solu- 
tions are found at the Alfven point where it takes the form 
of the Alfven regularity condition and at the base of the 
flow. This constrains the transfield equilibrium in that the 
Alfven regularity condition is imposed and the regularity 
of the magnetic surfaces at the Alfven critical surface is 
ensured. In our model, the outflow is parameterized by 
^{a), Q{a) and ao which is related to the mass loss rate. 
We deduce three flrst integrals of the motion, the speciflc 
angular momentum, the specific energy and the mass flux 
to magnetic flux ratio on a magnetic surface from criti- 
cality and Alfven regularity conditions. Different proflles 
for the variation of the entropy and the rotation rate with 
respect to the magnetic flux have been considered. The 
simpliflcations of the model allow to flnd analytical be- 
haviors of the flrst integrals as well as the shape of the 
critical surfaces in limiting cases. We found good agree- 
ment between analytical and numerical solutions. For a 
given entropy, magnetic rotators can be characterized by 
the ratio of the rotation rate and the magnetic to mass 
flux ratio on the polar axis, i.e. fl/ao- This latter ratio is 
given by the boundary conditions. 

Given the properties of the central emitting object, 
the model allows us to compute its corresponding dimen- 
sionless rotation parameter uj. Rotators can be deflned as 
slow, intermediate or fast according to whether ui is much 
less or close to unity or near its maximum value for fast 
rotators, (|)2. 

Slow rotators have quasi-spherical critical surfaces and 
a fast surface close to the Alfven surface, their properties 
strongly depend on the heating. 

Fast rotators have distorted critical surfaces and a fast 
surface far from the Alfven surface (approaching infln- 
ity only in the zero temperature limit). Their properties 
strongly depend on the magnetic flux and the rotation rate 
yet always retain a slow rotator behavior near the axis of 
rotation. They display a slow mode that acquires the char- 
acter of a sound wave as the rotation rate increases and 
are limited, in the present conical model, by a; = | ^ . 

For all types of rotators, it is found using conical 
shapes that the mass loss rate has a lower limit for a 
given rotation rate. The strongest effects on solutions are 
due to the rotation with respect the thermal and mass 
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loss rate effects for typical astrophysical values. Given the 
angular velocity fl{a) and the specific entropy Q(a), the 
solutions for the last first integrals of the motion, namely 
the specific energy E, the specific angular momentum L 
and the mass to magnetic flux ratio a can be determined 
numerically and in some limiting cases analytically. Non- 
constant rotation cases allow two- velocities outflow types, 
with possibly more complex solutions for different profiles 
of angular velocity and entropy. 

This simplified model makes it possible to investigate 
the structure of outflows far from the magnetized rota- 
tor source, unconstrained by type of boundary conditions, 
and without the need for self-similar assumptions. How 
outflows extend into the asymptotic region constitutes the 
subject of the stud y of P aper II (Lery et al.( |l998b| ). Pa- 
per III (Lery et al.( 1998c) will deal with a linear stability 
analysis of the asymptotic equilibria. 

Appendix A: Cold rotators 

In the conical case it useful to have the simplifled equa- 
tions for /3 <C 1 since numerically we have seen that this 
assumption is quite acceptable as long as Q(a) is not to 
big. These equations allow easy analytical works. Particu- 
larly for fast and slow rotators, one just have to study the 
limit of variation of uj to get the equations for positions 
and densities for the critical points. Thus we have 



Vs = 



Vf 



/3(7-l)x; 



2 1 + ^22.2(1 




(4-))^ 



(A.l) 



(A.2) 



One can deduce a reduced form for the relative energy and 
the relative gravity that only depend on Xs , Xf and on two 
parameters w, /3 
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Appendix B: Profiles of Q,{a) and Q{a) 

Definitions of parameters: 

• D,^: angular velocity 

• a: magnetic fiux 

• TV: large integer 

• B and C: constants 

• Ci and C2 ■ constants with Ci + C2 < 1 

Differential rotators 

— Rigid body rotation 

— Keplerian rotation 



(A.6) 



(B.l) 



Q = no{ 



(l-a^a(2-a))3/2 
l + exp{~Na + N/10) 

1 . 

(l + cxp(-7Va + iV/10))' 



(B.2) 
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Stellar type rotation 

n^floil-B sin^ e-C sin^ 6) 

with 

9 — arccos ( y/a{2 — a) J 

Decreasing stellar type rotation 
fl^floil-B cos^ 9-C cos'^ 9) 
Two velocities rotation (r^o-O) 

no 



fl = 



exp(iV(a-i)) + l-exp(-f) 
Two velocities rotation (fJo-^i) 

^ C2 



Q = no\Ci 



exp{N{a- i)) + 1 -exp(- 



N^ 



(B.3) 
(B.4) 

(B.5) 
(B.6) 

(B.7) 



Variations of Q{a) 
— Constant entropy 



= '^0 



(B. 



— Two entropies (Qo — ^' 



Q = Oo 0.5 



0.5 



N. I (B-9) 



exp(iV(a-i)) + l-exp(-f)^ 

— Increasing entropy 

Q = Qo(l-Ssin2 6'-Csin'*6') (B.IO) 

with 

9 = arccos (^/{2-a)) (B.ll) 

— Decreasing entropy 

n ^ rioil - B cos^ 9 - C cos'^ 9) (B.12) 

Appendix C: Examples of input parameters for 
astrophysical objects 
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